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Abstract 

Some recent developments on integrated damping mechanics for unidirectional 
composites, laminates, and composite structures are reviewed. Simplified damping 
micromechanics relate the damping of on-axis and off-axis composites to constituent 
properties, fiber volume ratio, fiber orientation, temperature, and moisture. Laminate 
and structural damping mechanics for thin composites are summarized. Discrete layer 
damping mechanics for thick laminates, including the effects of interlaminar shear 
damping, are developed and semi-analytical predictions of modal damping in thick 
simply-supported specialty composite plates are presented. Applications demonstrate the 
advantages of the unified mechanics, and illustrate the effect of fiber volume ratio, fiber 
orientation, structural geometry, and temperature on the damping. Additional damping 
predictions for composite plates of various laminations, aspect ratios, fiber content, and 
temperature illustrate the merits and ranges of applicability of each theory (thin or thick 
laminates). 


‘NASA Resident Research Associate at Lewis Research Center. 


INTRODUCTION 


The significance of passive damping to the dynamic performance of structures is widely 
recognized. Damping is a significant dynamic parameter for vibration and sound control, 
fatigue endurance, and impact resistance. Current light-weight and high-performance re- 
quirements imposed on most structural applications restrict the use of many traditional 
non-structural damping sources. Polymer- matrix composites are known to exhibit signif- 
icantly higher material damping than most common metals, as a result of the polymer 
matrix and their heterogeneity. Composites are already preferred in many structural ap- 
plications due to their high specific stiffness and strength, hence, the option of passive 
structural damping is an added advantage to these materials. 

Various damping mechanics theories for unidirectional composites and laminates [1- 
14] have been reported. Work on the damping of composite beams and plates [15-16], and 
plate/shell composite structures of general geometry and lamination has been also reported 
[17]. The research has demonstrated that composite damping depends on an array of 
micromechanical and laminate parameters, temperature, moisture, and existing damage. 
The composite damping is also anisotropic, but exhibits an opposite anisotropy trend 
than stiffness and strength, being minimum in the direction of the fibers and maximum 
in the transverse and shear directions. In addition, the structural composite damping 
was found to be strongly dependent on structural configuration and deformation state. 
It is apparent, therefore, that in order to realize significant structural benefits from the 
damping of composite materials, integrated damping mechanics are required, correlating 
the damping of composite structures to parameters of the basic constituent materials, 
laminate configuration, hygro-thermal conditions, and structural geometry. This paper 
reviews the past and present research at the Lewis Research Center in the development of 
integrated damping theories for thin and/or thick composites and composite structures. 
The damping mechanics complement previous work performed at Lewis on integrated 
composite mechanics [18]. 

These integrated methodologies enable the synthesis of global damping capacity at the 
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structural level, or any intermediate material level, while including the effects from crit- 
ical parameters at each individual material level (micromechanical, laminate, and struc- 
tural), friction from broken/debonded fibers, interlaminar layers, temperature, moisture, 
and global structural geometry and boundary conditions. 

Damping micromechanics explicitly relate the on-axis and off-axis damping capacities 
of the unidirectional composite with the micromechanical parameters and the fiber orienta- 
tion angle. Two different laminate/structural damping theories have been developed. The 
first was based on the classical laminate plate theory (CLPT) Kirchoff-Love’s assumptions 
which assume a linear displacement field through the thickness and neglect interlaminar 
shear effects. The CLPT damping theory combines simplicity, computational speed, and 
has yielded excellent results for thin laminates and laminated structures subjected predom- 
inantly to extensional and/or flexural deformations. An important element of the CLPT 
laminate damping theory is the inclusion of the interlaminar matrix layer damping, which 
proved to be a potentially important source of damping for angle-ply laminates. A finite- 
element based methodology for the synthesis of the passive damping and the simulation of 
the damped dynamic response of thin composite structures has also been developed. 

Recently, a novel discrete laminate damping theory for the damping of thick com- 
posites is developed. The method assumes a piecewise continuous displacement field with 
variable degrees of freedom through the laminate thickness. The effects of interlaminar 
shear damping are also included. The discrete damping theory, although computation- 
ally more cumbersome, is particularly suitable for thick composite laminates, or laminates 
where the representation of local interlaminar shear effects are important. A semi- analyt- 
ical methodology for the prediction of the dynamic properties of simply- supported (SS) 
specialty thick composite plates based on the aforementioned theory is presented, and re- 
sults for the modal damping values are reported. Comparisons with structural damping 
predictions based on the CLPT damping theory illustrate the advantages of both laminate 
theories. 

Application studies and experimental correlations illustrate the merits and the accu- 
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racy of these unified damping theories. In addition, they provide valuable information and 
guidelines for the design of composite laminates/structures with good dynamic character- 
istics, as well as, for the successful measurement of damping in composites. 

COMPOSITE DAMPING 

This section briefly reviews the synthesis of damping for on-axis (damping along the 
material axes) and off-axis composites (unidirectional composites loaded at an angle), 
including the effects of temperature and moisture on composite damping. Additional 
details are provided in Ref. [13] . 

On- Axis Damping 

For a unidirectional composite loaded along the material axes (Fig. la), closed-form 
expressions have been developed for the synthesis of damping based on elastic hysteretic 
damping assumptions [13]. For orthotropic but transversely isotropic materials, like uni- 
directional reinforced composites, six damping coefficients will completely characterize the 
damping of the composite, that is: (1) longitudinal damping V’/i (direction 11), (2) trans- 
verse in- plane damping (direction 22), (3) transverse through-the- thickness damping 

(direction 33), (4) in-plane shear damping (direction 12), (5) interlaminar shear 
damping V’M (direction 23), and (6) interlaminar shear damping V’is (direction 13). These 
six damping capacities are explicitly related to fiber/matrix moduli and damping values, 
and the fiber volume ratio (FVR). Typical damping predictions and experimental correla- 
tions for Gr/Epoxy unidirectional composites are shown in Fig. 2. 

Off-Axis Damping 

For the case of off-axis composites, ie. composites loaded at an angle (Fig. lb), a 
transformation provides the equivalent damping capacity of the composite in the structural 
coordinate system. Ref. [13] presents the damping transformation for in-plane SDCs. The 
damping transformation is extended herein, based on the same principle of invariance of 
the dissipated energy, to include out-of-plane shear and transverse damping. In such case, 
the off-axis composite damping is best described by the following 6 by 6 damping matrix 
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[V'c] = [*] T [M*]" T 


( 1 ) 


[•0c] 


The transformation matrices [il] T and [R]~ T are shown in the Appendix and induce the 
effect of fiber orientation angle, while the diagonal matrix [V»i] represents the on-axis damp- 
ing values 
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( 3 ) 


The non-diagonal terms 12 indicate coupling between in-plane extensions, the terms 16 and 
26 coupling between the in-plane axial and shear stresses, while the terms 45 the coupling 
between interlaminar shear stresses. Apparently, off-axis loading will affect the overall 
damping capacity of the composite in two distinct, yet uncommon to isotropic materials, 
ways: (1) by altering the values of the diagonal terms, which is equivalent to altering the 
dissipative capability of the ply directly associated to normal and shear strains, and (2) 
by inducing and altering non-diagonal terms, which control the amount of strain energy 
dissipated by coupled deformation modes. This is best illustrated in Fig. 3 where the 
equivalent axial damping, and the values of the damping terms V’cii? ^22? and ^66 of a 
0.50 FVR Gr/Epoxy composite vs. the fiber orientation are plotted. The important effects 
of fiber orientation and coupling are apparent. 
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Hygrothermal Effect 

Many mechanical properties of polymers, including damping, are sensitive to hy- 
grothermal variations. In this manner, temperature and moisture may have a definite 
effect on the response of polymer composites. Previous studies [18,19] have shown that 
the hygro-thermal effect on most mechanical properties of the matrix can be modelled as, 

Pm r Tgw — T_ i o .5 /•4I 

Po [ T gd -Tj {) 

where: P and T represent matrix properties and temperature; subscript M represents 
Mechanical Properties; and subscripts o, gd, gw respectively indicate room, dry glass 
transition and *wet, gluss transition conditions. An inverted form of eq. (4) has been 
proposea for the hygro-thermal effect on the matrix damping: 


Pi 

Po 



The wet glass transition temperature is [19], 


( 5 ) 


Tgw = TgdiO.OOomt - 0.1m/ + 1) (6) 

where m/ is the moisture content. The exponent q in eq. (5) can be correlated to experi- 
mental damping data of each individual polymer matrix. In the present paper an average 
value q=0.5 is assumed. The effects of temperature and moisture on the composite damp- 
ing are introduced based on the previously described micromechanics. A similar approach 
is used for other mechanical properties. 


LAMINATE DAMPING 

Two different laminate damping mechanics theories have been developed for thin and 
thick laminates, respectively. The first is based on the CLPT Kirchoff-Love assumptions 
and has provided excellent correlations with experimental results. The second is based 
on discrete layer laminate plate theory (DLPT) and incorporates a piecewise continuous 
deflection field of variable degrees of freedom. The discrete layer damping theory (DLDT) 
is developed to model the damping of thick laminates, or laminates where the interlaminar 
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shear effects are non-negligible. The assumed displacement fields for both laminate theories 
are shown in Fig. 4. In both cases, the dissipated strain energy in the laminate is defined 
as, 

r k/2 

Aw l = 1/2 / e c T [E c ][rJ> c }e c dz (7) 

J-h / 2 

and the maximum laminate strain energy 

fh/2 
-h/2 

where [E c ] is the off-axis stiffness matrix (see Appendix). The equivalent laminate specific 
damping capacity (SDC) then is, 


r h/2 rr 

wl = 1/2 / e c T [E c ]e c dz 
J-h! 2 


( 8 ) 


= Awl/w l 


( 9 ) 


Thin Laminates 

An overview of the laminate damping theory for thin laminates [14] is provided here. 
Based on the Kirchoff-Love’s assumptions of uniform linear displacements through the 
thickness and no interlaminar shear, the laminate strain e c contains only in-plane strains 
which are related to the mid- plane strain e° and curvature k as follows: 

e ci = {e°,k l }{l,z} T (* = 1,2,6) (10) 


Combination of eqs. (7,10) and integration through the thickness provides the dissipated 
strain energy of the laminate, expressed in terms of the 3 by 3 damping matrices [Ad], [Cj], 
and [Dd] (see Appendix) which are related to the properties of the composite laminate only. 


Aw l = ^{ e °> k > 


[^<] [cy 

[Ci] [Dd] 



( 11 ) 


The previous expression is complete, in that, it represents general laminate configurations 
and deformations. Indeed, the extensional damping matrix [A^] represents the laminate 
damping capacity in pure extension, the flexural damping matrix [Dd] represent the lam- 
inate damping in pure bending and torsion, and the coupling damping matrix [Cj] is 
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associated with coupled extension-flexure deformations. All three matrices in eq. (11) are 
fully populated. They also include the interlaminar damping contributions of the interply 
matrix layer, which have shown to be a significant additional source of damping in compos- 
ite laminates [14]. Similarly, the maximum laminate strain energy in eq. (8) is expressed 
as: 


w -LoufW w]i e °\ 

WL - 2 {e ,k} |j£] 


( 12 ) 


where [A], [C], and [Z>] are the extensional, coupling, and flexural laminate stiffness ma- 
trices. Therefore, the equivalent damping capacity of a general laminate is calculated by 
combining eqs. (9,11,12) and is a function of the composite laminate parameters and the 
specific deformation state {e°,k}. 

This laminate damping theory has yielded excellent correlations with experimental 
results. Typical damping predictions for a 0.50 FVR graphite/epoxy (0 4- 0/90 + 0/45 + 
0/ — 45 + 0), laminate under cyclic free-flexure (a bending moment M = (M x ,0,0)) is 
shown in Fig. 5. Additional laminate damping predictions and correlations for various 
laminate configurations are shown in ref. [14]. 


Thick Laminates 


A new discrete-layer laminate damping theory (DLDT) is recently developed for thick 
laminates and/or laminates with significant interlaminar shear effects. In contrast to the 
CLPT, the present theory assumes a general piecewise continuous in-plane displacement 
field through the laminate thickness, in addition to the midplane deflections. Discrete 
layer theories were historically proposed by Grigolyuk [20] and subsequently generalized 
(see review [21]) for the more accurate calculation of stresses in thick laminates. The 
authors, however, have realized that they combine the potential for accurate damping 
predictions in thick laminates, or laminates with weak interlaminar damping layers, where 
deviations from the CLPT assumptions are expected. The assumed displacement field has 
the form, 

u(x,y,z,t) = u°(x,y,t) + t i(x,y,z,t) 
v(x,y,z,t) = v°(x,y,t) + v(x,y,z,t) 
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(13.1) 

(13.2) 


w{x,y,z,t) = w°(x,y,t) (13.3) 

where superscript o represents the uniform through-the-thickness midplane deflection. Ap- 
plying a piecewise continuous linear displacement field, as suggested by Reddy and co- 
workers [22-23], the previous eqs. take the form: 

N 

u (x,y,z,t) = u°(x,y,t ) + ^11’ (x,y,i)F 3 (z) (14.1) 

3 = 1 

N 

v(x,y,z,t) = v°(x,y,t) + v j (x, y, i) F j (z) (14.2) 

3=1 

w(x,y,z,t) = w°(x,y,t) (14.3) 

where, u J and u* 7 are displacements, along the x and y directions respectively, preferably 
at the interfaces between various plies, or sublaminates (group of plies), and F J (z) are 
linear interpolation functions. In this manner, the assumed in-plane displacement field is 
general, in that it may represent general displacement fields and interlaminar shear strains 
through-the-laminate thickness in the laminate plane. 1 

The laminate strains are directly derived from eqs. (14) 

N 

e a = e °ci + '52 e ii F3 ( z ) i = 1,2,6 (15.1) 

3 = 1 
N 

£ « = t ? i = 4 ’ 5 < 15 - 2 ) 
3=1 

where, the midplane strains are, 
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and the generalized strains are, 


(16.1) 
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(17.1) 
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4 = v ’ 4 = ui (17.2) 

The comma in the subscripts indicates differentiation. Combination of eqs. (7, 15-17) 
ultimately provides the dissipated strain energy per unit area, 

Aw l = l/2(«f \AiV' + 2 ef £>>’ + '£'£ efloriC) (18) 

j=l j=l m=l 

where, [Ad] is the previously defined extensional laminate damping matrix, but with ad- 
ditional out-of-plane shear terms. The generalized coupling damping matrices [B d ] and 
flexural/ shear matrices [D 3 ™] are new. 


Ni *hk 

\A d ] = E / {EM'l’Mz 


(S’) in = (IJE'UMM*. F’{z)d. 

k=l Jh k-x 
k=l 

(cr)in = E r ({EMi’M™F i U)F m (z)d 

i=l Jh k-1 

(Din in = E r ((EM'l’'}k)<«F!,(z)F;(z)i 

k=l Jh k-1 


z i,n = 1,2,6 
z t,n = 4,5 
z i,n = 1,2,6 


(19.1) 


(19.2) 


(19.3) 


z i,n = 4, 5 


In a similar manner, the combination of eqs. (8. 15-17) provides the maximum strain 
energy per unit area, 

N N N 

w L = l/2(*f {Atf + 2<f El^K + E E 'f i^”)C) (20) 

j=l j=l m=l 

where, [^4] is the previously defined extension laminate damping matrix, however, with 
additional out-of-plane shear terms. The generalized coupling damping matrices [B J ] and 
flexural/shear matrices [D^ m ] are new. 
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i,n = 1,2,6 


( 21 . 2 ) 


(B’U = E r (\E'] k )i„F>(z)dz 

k=l Jh k-1 


(B>) 


N ‘ rh 
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k=l Jh *-' 
* rh k 
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(D im )<„ = E / ‘ ([^]t)inF i (^)F m (^)<i 
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(21.3) 


t, n = 4, 5 


The equivalent laminate damping for a given local displacement field is provided by 
eq. (9). 


STRUCTURAL DAMPING 


The present section describes the procedure followed for characterizing the global 
damping capacity of composite structures. It is clear from the previous sections, that the 
damping capacity of the composite structure will strongly depend on the global deforma- 
tion state. This difficulty is overcome by considering the modal damping values, that is, 
the damping of the structure vibrating solely at the respective modal shapes. Since the 
vibration modes form a complete set representing the elasto-dvnamic response of the struc- 
ture undergoing small deflections, the modal SDCs provide a complete characterization of 
the global structural damping. The modal damping coefficients, in connection with other 
modal parameters, can also provide the damped dynamic response of the structure [24]. 

Thus, the objective of this section is to present the method for synthesizing modal 
damping. The modal SDC associated with the n-th vibration mode Is: 


, _ JyAw n dV 

^ JyW n dV 


( 22 ) 


where Au; n and w n represent the dissipated and maximum specific modal strain energies 
per cycle respectively. For the case of laminated composite structures the structural modal 
damping becomes, 


n = 


f A ^WLndA 
f A WL n dA 


(23) 
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where A wi n and wi n are respectively the dissipated and maximum laminate modal strain 
energies of the n-th mode. 

Thin Composite Structures 

For thin composite structures, eq. (23) is combined with eqs. (11-12) to yield the 
modal damping. For specialty laminates and simple structural configurations, such as 
plates and shells, exact mode shapes may be found, and explicit integration of eq. (23) 
is possible, resulting in analytical expressions for modal damping. Such expressions have 
been derived by the authors for simply-supported rectangular plates. 

More general laminate/structural configurations, however, require a form of approxir 
mate solution. As a result, approximate procedures based on finite element discretization 
were developed [17]. The integrations in eq. (23) are performed numerically, first over 
the area of the element, and then over the structural area. A specialty finite element was 
developed with 6 degrees of freedom per node. The damping capacity of the element was 
represented by the element damping matrix and included membrane, flexural, and coupling 
damping terms. The complete inclusion of extensional, flexural and coupling terms into 
the damping element matrix, readily provided the capacity for predicting the damping in 
composite structures with general laminate configurations and/or material and geomet- 
ric coupling (for example, thin composite shells, thin- wall box or tubular beams, and so 
forth). Typical predicted SDCs for three graphite/epoxy composite structures of varying 
structural complexity, a cantilever beam, plate and shell respectively, are shown in Fig. 6. 
The assumed lamination was (±0)s,. Among other issues, Fig. 6 dramatically illustrates 
the effect of fiber orientation, structural aspect ratio, curvature, and mode order on the 
structural modal SDCs. 

Thick Composite Plates 

A new semi-analytical method for the calculation of damping in simply- supported 
thick cross-ply composite plates, in connection with the discrete layer laminate damping 
theory previously described. Assuming a rectangular a by (3 symmetric composite plate 
with negligible coupling (.4i6=A26=0, B{ 6 =b1 6 =0, D{™ =D^™= 0), the following Navier 
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fundamental solutions form a complete set of mode shapes in the x-y plane: 



(24.1) 

= '^ n 5m(ai)co5(i>y)e’ a ’ t 

(24.2) 

w mn( x ^y^) = W^ n sin(ax)sin(by)e xu,t 

(24.3) 


(24.4) 

= l / mn sin ( aa: ) coa ( 6 y) eta ’ 1 

(24.5) 


Where, a = m7r/a and b = n7r//3. Combination of eqs. (16,17,24) yields the modal 
mid-plane and generalized strains as separable functions of x, y coordinates, time, and 
amplitudes. 

{ c „> = [*•.]{ = ibu^L}^ ( 25 ) 

where, the terms in matrices [i?] mn are sine/cosine functions of x, y coordinates and mode 
order. The amplitude displacement vectors are U° = {£7°, V° , W°} T and U J = {C/ J , V*} T . 
The kinetic energy through the thickness of the plate is: 

rh/2 

-h/2 

Considering eqs. (14) the laminate kinetic energy per unit area takes the form: 

k l = i/2({i»} r [4 M ){u «} + + e e m 

j= 1 j = 1 m= 1 

where the generalized laminate mass matrices are 

r^k 

[Am ) = diag(p k )dz 

k= 1 


fh/2 

Kl = 1/2 / {u} T p{u}dz 

J-b.ii 


(26) 


(28.1) 


\B’ m ] = E / diag(f> t )F’ (z)iz 
k= 1 


( 28 . 2 ) 
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diag(p k )F 3 (z)F m (z)dz 


(28.3) 



The term diag(pk) indicates a diagonal matrix, with all diagonal terms equal to the 
density of the k-th ply. 

By substituting eqs. (25) into the strain and kinetic laminate energy, eqs. (20) and 
(27), integrating over the plate area, and applying Lagrangian dynamics, it is proved that 
the undamped modal analysis (free vibration) solution of the plate takes the form: 

-u4 n [Jl/ mn ]U mn + [A' mn ]u mn = 0 (29) 

■where, U mr - {U^ n ; U^ n , ...,U^ n }. Numerical solution of this eigenvalue problem 
provides the natural frequencies and the through-the-thickness modes for each order mn 
of plane modes, in the context of eqs. (24). 

The modal damping of each mode is subsequently calculated, in accordance with eq. 
(23). First eqs. (19) and (21) are integrated through-the-thickness, and are combined with 
eqs. (18) and (20) respectively. Then the generalized stresses and strains are related to 
the calculated mode shape using eq. (25). The integrations over the area of the plate are 
performed analytically. Only the natural frequencies and the mode shapes through-the- 
thickness are calculated numerically, hence, the whole procedure is semi- analytical. 

An analogous procedure was also developed for the exact analytical solution of the 
natural frequencies and modal SDCs of thin rectangular plates based on the CLPT as- 
sumptions. In the remaining paragraphs, predictions wfith both thick and thin composite 
damping theories are shown. 

Predicted results for the first 10 modes of a square 406 mm by 406 (16in by 16in) 
simply-supported (SS) plate, 4.06 mm thick (0.16 in), of 16 0.50 FVR Gr/Epoxy plies, 
each 0.01 in thick, based on both discrete layer plate theory (DLPT) and CLPT are shown 
in Table 1. Three laminate configurations are investigated: (1) a unidirectional plate (0)i6, 
(2) cross-plies ( 04 / 904 ),, and (3) alternating cross-plies (0/90)4,. The fibers in the 0 degs 
plies are parallel with the x axis and 16 discrete linear displacement segments, one in each 
ply, were used (N=17). As seen, the modal SDCs vary widely with the modal shape. 
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In all cases, the differences in the prediction of the first few modal SDCs is small and 
both theories show good agreement in the prediction of natural frequencies. For this plate 
of high aspect ratio, this was expected and lends credence to the method. However, the 
differences between the modal damping of higher order modes of the 0/90 plates, is greater 
than in the case of the unidirectional plate. This is attributed to the higher interlaminar 
shear stresses developed in the transverse plies in order to balance the normal stresses in 
the longitudinal plies. 

The fundamental modal SDC and natural frequency of a square 0.50 FVR (04/904), 
Gr/Epoxy plate are plotted in Fig. 7 for varying aspect ratios a/h. At aspect ratios 
greater that 50, both CLPT and DLPT give comparable predictions, which again rein- 
forces the validity of the DLPT damping method. At aspect ratios lower that 50, the 
DLPT damping method gives significantly higher predictions of damping than the CLPT, 
because the interlaminar shear effects are important, which illustrates the superiority of 
the methodology. The differences in frequency predictions are small. 

To illustrate the potential of the unified composite damping methodologies presented 
herein, the fundamental SDC and natural frequency of a 203 mm by 203 mm (8 in by 8 
in), (04/904), Gr/epoxy plate for varying FVR, either in the 90 degs plies only or the 0 
degs plies only, are shown in Fig. 8. The modal damping is rather insensitive to the FVR 
variation in the 90 degs plies, however, the DLPT predicts significantly higher damping at 
low FVR. It is recalled that at low FVR, the shear modulus G 123 is reduced, hence, higher 
strain energy is stored in the 90 degs plies to balance the bending stress in the 0 degs plies, 
moreover, the interlaminar shear damping is increased. The variation of the FVR in the 
0 degs plies has a definite effect on the modal damping and frequency, but both theories 
provide equivalent results. 

Predictions of the fundamental SDC and natural frequency of the same plate, sub- 
jected to uniform through-the-thickness temperature variations are plotted in Fig. 9. The 
CLPT significantly underestimates the fundamental SDC at higher temperatures. It is 
recalled, that in higher temperatures, the shear modulus of the 90 degs plies decreases and 
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the shear damping increases respectively. The temperature variation has little effect on 
the natural frequency. 

CONCLUSIONS 

The paper presented an update on past and present research work performed by the 
authors on the development of unified theories for the damping of polymer composites and 
composite structures at NASA-Lewis Research Center. Damping mechanics for compos- 
ites were developed in a unified sense, that is, for all three stages of structural complexity 
(micromechanical, laminate, and structural), including hygrothermal effects, and effects of 
interlaminar shear. Laminate and structural damping mechanics for thin composite lami- 
nates of general laminate configurations were first developed and reported. The damping 
mechanics can handle exteusional, flexural and coupled deformations, therefore, can pre- 
dict the damping of symmetric, antisymmetric, and/or fully asymmetric laminates, or the 
damping of co^' ‘.e structures with geometric coupling such as thin-wall beams, shells, 
and sc lorfch. Representative results and experimental correlations were presented. 

Damping mechanics for thick composite laminates/structures are presently under de- 
velopment and a unified discrete-layer damping theory was presented. The discrete dis- 
placement field of variable degrees of freedom enabled the calculation of interlaminar shear 
strains and the incorporation of interlaminar shear damping. Semi-analytical solutions 
were developed for the prediction of damped dynamic characteristics of thick composite 
laminates. Comparisons with the CLPT damping theory illustrated that for high aspect 
ratios both theories give comparable results. However the discrete- layer theory provided 
better predictions, in specific cross-ply laminate configurations, low aspect ratios, high- 
order modes, cross-plies with inner plies of low FVR, and cross-plies at higher tempera- 
tures. Both theories yielded comparable predictions of natural frequencies in most cases. 

Overall, the presented evaluations of the method demonstrated the power and ver- 
satility of these integrated composite damping mechanics. The damping mechanics have 
been incorporated into in-house research codes, including ICAN [ 18 ]. Subjects of future 
research will be composite laminates with interlaminar damping layers, and the effects 
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of damping improvements on the stiffness, natural frequencies, buckling resistance, and 
impact performance of thick composite structures. 
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APPENDIX 


Off-axis ply stiffness matrix {E c } (ref. 18): 


{crj = [E c ]{eJ 


(Al) 


[E c ] = [R]- 1 [E,][R]- t 


(A2) 


Transformation matrices: 


' 771 2 71 2 0 0 0 2 77171 “ 

71 2 771 2 0 0 0 —277171 


[*]= 0 0 0 m! 0 ’ n = sine ( A 3) 


0 0 0 — n m 0 

,-mn mn 0 0 0 m 2 — n 2 . 


[R]- 1 = [*(-*)] 


(A4) 
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Laminate damping matrices (CLPT): 


N , Ni- 1 

[Acd] = Y M&lfhMi + Y 
1=1 1=1 

[c CD ]=x; t( 2 f - + Y 

1=1 i= 1 

i£ CD ] = f; i( 2 ?-4) i [£;c]n«+ y 4. i-ffiiw-ii 


(> 15 ) 

(A6) 

(A7) 


t=l t=l 

The interply distortion energy coefficient H, and matrix [S] are described in ref. 18 . 
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Table 1. Modal characteristics of square simply supported plates 


Mode ( 11 , 111 ) 


Nat. frequency (Hz) 
(DLPT) (CLPT) 


Modal Damping, % 
(DLPT) (CLPT) 


(<W Plate 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


1 

2 

3 

4 


10 


(1.1) 

128.8 

129.2 

1.035 

0.9% 

(1.2) 

161.5 

161.9 

2.561 

2.531 

(U) 

240.2 

240.8 

3.861 

3.842 

(1,4) 

3673 

368.5 

4.399 

4.389 

(2.1) 

490.8 

496.8 

0.761 

0.606 

(2,2) 

510.6 

516.8 

1.148 

0.9% 

(1,5) 

538.8 

541.2 

4.582 

4.575 

(23) 

556.3 

562.8 

1.827 

1.684 

(2,4) 

640.7 

647.8 

2.654 

2.531 

(2,5) 

772.3 

780.4 

3.390 

3.294 

s Plate 

(i,i) 

128.8 

129.2 

1.026 

0.9% 

(1.2) 

231.1 

231.9 

1.464 

1.428 

(13) 

449.5 

452.8 

1.440 

1.357 

(2,1) 

461.7 

468.2 

0.764 

0.635 

(2,2) 

510.1 

516.8 

1.113 

0.9% 

(23) 

652.5 

660.5 

1.437 

1.326 

(1,4) 

767.6 

777.5 

.. 1.428 

1.279 

(2,4) 

914.4 

927.7 

1.572 

1.428 

(3,1) 

1011.0 

1043.0 

0.845 

0.562 

(3,2) 

1041.0 

1074.0 

1.012 

0.736 

Plate 

(i,i) 

128.9 

129.2 

1.019 

0.9% 

(1,2) 

338.0 

340.3 

0.931 

0.865 

(2,1) 

392.9 

396.5 

0.823 

0.737 

(2,2) 

512.1 

5.168 

1.076 

0.9% 

(13) 

720.8 

732.0 

0.904 

0.752 

(23) 

827.2 

839.6 

1.103 

0.964 

(3,1) 

852.1 

869.7 

0.837 

0.643 

(3,2) 

922.2 

940.3 

1.025 

0.846 

(33) 

1140.0 

1163.0 

1.173 

0.9% 

(1,4) 

1254.0 

1288.0 

0.972 

0.706 
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FIGURE 1 —UNIDIRECTIONAL COMPOSITES. 
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FIGURE 2 — TYPICAL SDCS FOR HM-S GRAPHITE/POLYESTER. 
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FIGURE 4 — KINEMATIC ASSUMPTIONS IN BOTH COMPOSITE THEORIES. 
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FIGURE 3 — EFFECT OF RBER ORIENTATION ON 
THE DAMPING OF A 0.50 FVR KM-S/EPOXY 
COMPOSITES 



FIGURE 5 — FLEXURAL DAMPING OF A 0.50 FVR HM-S/EPOXY 
(ft«)+a/45+8/-45*6) s LAMINATE 
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FIGURE 6 -EFFECT OF PLY ANGLE 6 ON THE MODAL PROPERTIES OF VARIOUS SYMMETRIC ANGLE-PLY GRAPHITE/EPOXY COMPOSITE STRUCTURES. 
DIMENSIONS ARE IN INCHES 
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ASPECT RATIO (a/ti) 


FIGURE 7 — EFFECT OF ASPECT RATIO ON THE FUNDA 
MENTAL MOOAL CHARACTERISTICS OF A (Q, /9C^ )s 
HM-S/EPOXY PLATE. 
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FIGURE 8 — EFFECT OF FVR ON THE FUNDAMENTAL MOOAL 
CHARACTERISTICS OF AN (0 4 /9Q<)s HM-S/EPOXY PLATE 
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FIGURE 9 -EFFECT OF TEMPERATURE ON THE FUNDA- 
MENTAL MOOAL PROPERTIES OF A (0 4 *0 4 )s 
HM-S/EPOXY PLATE 
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